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Abstract 

We calculate the in-medium masses of J/tp and of the excited states of charmonium (^(3686) and 
7/^(3770)) in isospin asymmetric nuclear matter at finite temperatures. These mass modifications 
arise due to the interaction of the charmonium states with the gluon condensates of QCD, simulated 
by a scalar dilaton field introduced to incorporate the broken scale invariance of QCD within an 
effective chiral model. The change in the mass of J/ip in the nuclear matter with density is seen 
to be rather small, as has been shown in the literature by using various approaches, whereas, the 
masses of the excited states of charmonium (i/;(3686) and '0(3770)) are seen to have considerable 
drop at high densities. The dependence of the masses of the charmonium states on the isospin 
asymmetry has also been investigated in the hot nuclear matter and is seen to be appreciable at 
moderate temperatures and high densities. These medium modifications of the charmonium states 
should modify the experimental observables arising from the compressed baryonic matter produced 
in asymmetric heavy ion collision experiments in the future facility of FAIR, GSI. 
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I. INTRODUCTION 



The study of the in-medium properties of hadrons is an active topic of research in strong 
interaction physics, both experimentally and theoretically. The topic is of direct relevance 
in the context of heavy ion collision experiments, which probe matter at high temperatures 
and/or densities. The medium modifications of the hadrons have direct consequences on 
the experimental observables from the strongly interacting matter produced in heavy ion 
collision experiments. The in-medium properties of the kaons and antikaons in the nuclear 
medium are of relevance in neutron star phenomenology where an attractive interaction of 
antikaon-nucleon can lead to antikaon condensation in the interior of the neutron stars. The 
medium modifications of kaons and antikaons also show in their production and propagation 
in the heavy ion collision experiments. The modifications of the properties of the charm 
mesons, D and D as well as the J / ip mesons and the excited states of charmonium, can have 
important consequences on the production of open charm and the suppression of J/ip in 
heavy ion collision experiments. In high energy heavy ion collision experiments at RHIC as 
well as LHC, the suppression of J ftp can arise from the formation of the quark-gluon-plasma 
(QGP) [l|,|2|. 

The D (D) mesons are made up of one heavy charm quark (antiquark) and one light (u 
or d) antiquark (quark). In the QCD sum rule calculations, the mass modifications of D 
(D) mesons in the nuclear medium arise due to the interactions of light antiquark (quark) 
present in the D{D) mesons with the light quark condensate. There is appreciable change 
in the light quark condensate in the nuclear medium and hence D (D) meson mass, due to 
its interaction with the light quark condensate, changes appreciably in the hadronic matter. 
On the other hand, the charmonium states are made up of a heavy charm quark and a 
charm antiquark. Within QCD sum rules, it is suggested that these heavy charmonium 
states interact with the nuclear medium through the gluon condensates. This is contrary to 
the interaction of the light vector mesons (p, u, 0), which interact with the nuclear medium 
through the quark condensates. This is because all the heavy quark condensates can be 
related to the gluon condensates via heavy-quark expansion js | . Also in the nuclear medium 
there are no valence charm quarks to leading order in density and any interaction with the 
medium is gluonic. The QCD sum rules has been used to study the medium modification of 
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D mesons j4| and light vector mesons The QCD sum rule approach jo] and leading order 

perturbative calculations [3] to study the medium modifications of charmonium, show that 

ii 

the mass of J/ip is reduced only slightly in the nuclear medium. In [Sj, the mass modification 
of charmonium has been studied using leading order QCD formula and the linear density 
approximation for the gluon condensate in the nuclear medium. This shows a small drop 
for the J/ip mass at the nuclear matter density, but there is seen to be significant shift in 
the masses of the excited states of charmonium {ip{3Q86) and '?/'(3770)). 

In the present work, we study the medium modification of the masses of J/ip and excited 
charmonium states ?/'(3686) and ?/'(3770) in the nuclear medium due to the interaction with 
the gluon condensates using the leading order QCD formula. The gluon condensate in 
the nuclear medium is calculated from the medium modification of a scalar dilaton field 
introduced within a chiral SU(3) model 9| through a scale symmetry breaking term in the 
Lagrangian density leading to the QCD trace anomaly. In the chiral SU(3) model, the 
gluon condensate is related to the fourth power of the dilaton field x ^^id the changes in 
the dilaton field with density are seen to be small. The model has been used successfully 



to study t 
matter in 



le medium modifications of kaons and antikaons in isospin asymmetric nuclear 



10| and in hyperonic matter in [llj]. The model has also been used to study the 



D mesons in asymmetric nuclear matter at zero temperature [12| and in the symmetric and 
asymmetric nuclear matter at finite temperatures in Ref. 13| and [l^. The vector mesons 
have also been studied within the model jisl . 16|. In the present investigation, we study 
the isospin dependence of the in-medium masses of charmonium obtained from the dilaton 
field, X calculated for the asymmetric nuclear matter at finite temperatures. This study will 
be of relevance for the experimental observables from high density matter produced in the 
asymmetric nuclear collisions at the future facility at GSI. 

The outline of the paper is as follows : In section II, we give a brief introduction of 
the chiral SU (3) model used to study the in-medium masses of charmonium in the present 
investigation. The medium modifications of the charmonium masses arise from the medium 
modification of a scalar dilaton field introduced in the hadronic model to incorporate broken 
scale invariance of QCD leading to QCD trace anomaly. In section III, we summarize the 
results obtained in the present investigation. 
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II. THE HADRONIC CHIRAL 5C/(3) x SU{3) MODEL 



We use an effective chiral SU{3) model for the present investigation 9|]. The model is 



3ased on the nonlinear realization of chiral symmetry 



lJl9 



and broken scale invariance 



9], [l5|, ll6|. This model has been used successfully to describe nuclear matter, finite nuclei, 
hypernuclei and neutron stars. The effective hadronic chiral Lagrangian density contains 
the following terms 

BW + ^vec + -^^0 + ^SB (l) 

W=X,Yy,A,u 

In Eq. ([T]), Ckin is kinetic energy term, Cbw is the baryon-meson interaction term in 
which the baryon-spin-0 meson interaction term generates the vacuum baryon masses. C^ec 
describes the dynamical mass generation of the vector mesons via couplings to the scalar 
mesons and contain additionally quartic self-interactions of the vector fields. Cq contains 
the meson-meson interaction terms inducing the spontaneous breaking of chiral symmerty 
as well as a scale invariance breaking logarthimic potential. Csb describes the explicit chiral 
symmetry breaking. 

To study the hadron properties at finite temperature and densities in the present inves- 
tigation, we use the mean field approximation, where all the meson fields are treated as 
classical fields. In this approximation, only the scalar and the vector fields contribute to the 
baryon-meson interaction, Cbw since for all the other mesons, the expectation values are 
zero. The interactions of the scalar mesons and vector mesons with the baryons are given 
as 

Csscal + Cbvcc = - E ^» + QuiilO^ + Qpil^P + 5'<^i7O0] ^i- (2) 

i 

The interaction of the vector mesons, of the scalar fields and the interaction corresponding 
to the explicitly symmetry breaking in the mean field approximation are given as 

r ^ f 2 2 , 22, 2 i2\ X 

^ Xo 
+ (?4(c^' + 6a;V + p^ + 20^), (3) 



Co = -hox' {a' + e + 6') + h {cr' + (' + 6' 
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5S 
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The effective mass of the baryon of species i is given as 



(4) 



(5) 



(6) 



The baryon-scalar meson interactions, as can be seen from equation (|6]), generate the baryon 
masses through the couphng of baryons to the non-strange a, strange ( scalar mesons and 
also to scalar-isovector meson 6. In analogy to the baryon-scalar meson coupling there 
exist two independent baryon-vector meson interaction terms corresponding to the F-type 
(antisymmetric) and D-type (symmetric) couplings. Here antisymmetric coupling is used 
because the universality principle j2o| and vector meson dominance model suggest small 
symmetric coupling. Additionally, we choose the parameters so as to decouple the 

strange vector field 0^ ~ s'y^s from the nucleon, corresponding to an ideal mixing between 
oj and mesons. A small deviation of the mixing angle from ideal mixing 2ll-l23| has not 
been taken into account in the present investigation. 

The concept of broken scale invariance leading to the trace anomaly in (massless) QCD, 
9f: = ^^f^G°' ^uG^'^"' , where G^^, is the gluon field strength tensor of QCD, is simulated in the 



effective Lagrangian at tree level [2^ through the introduction of the scale breaking terms 



C 



scalebreaking 



1 4i IX 

-^X In 



Xo 



X_ 

Xo 



(7) 



,det(X)o 

where /s = det(X), with X as the multiplet for the scalar mesons. These scale breaking 
terms, in the mean field approximation, are given by the last two terms of the Lagrangian 
density, given by equation (jlj). The effect of these logarithmic terms is to break the scale 



invariance, which leads to the trace of the energy momentum tensor as 



25| 



0', 



x|^-4£ = -(l-rf)x^ 



(8) 



Hence the scalar gluon condensate of QCD {{C^ ^uG^^"")) is simulated by a scalar dilaton 
field in the present hadronic model. 

The coupled equations of motion for the non-strange scalar field cr, strange scalar field 
scalar-isovector field 5 and dilaton field Xi ^"^^ derived from the Lagrangian density and are 



given as 



kQX^a - 4fci ((T^ + C' + 5') ^ - 2/C2 (a^ + 3^52) _ 2hx<yC 



2a 



dxr , x_ 

3 C Uo 



(9) 



(10) 



(11) 



koX cr' + e + 6']~k, a'-6')C + X 



4 / ,t^2 x2 

6 \\ (Tom 



XoJ 



2x 

Xo 



+ (4^4 - d)x^ 
1 



In the above, pi'^ are the scalar densities for the baryons, given as 

d^k m* ( 1 1 



Pi = li 



+ 



0(12) 



(13) 



(27r)3 E*{k) \ e(^'*W-/^'*)/^ + 1 (,{Ei'{k)+in*)/T j^i 

where, E*{k) = {k"^ + m*'^)^/'^^ and, p* = pi — g^itu — Qpip — g^fnCp, are the single particle 
energy and the effective chemical potential for the baryon of species i, and, 7^=2 is the spin 
degeneracy factor [lo| . 

The above coupled equations of motion are solved to obtain the density and temperature 
dependent values of the scalar fields (a, ( and 6) and the dilaton field, in fhe isospin 
asymmetric hot nuclear medium. As has been already mentioned, the value of the x is related 
to the scalar gluon condensate in the hot hadronic medium, and is used to compute the in- 
medium masses of charmonium states, in the present investigation. The isospin asymmetry 
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in the medium is introduced through the scalar-isovector field 6 and therefore the dilaton 
field obtained after solving the above equations is also dependent on the isospin asymmetry 
parameter, rj defined as r] = {pn — Pp)/{2pB), where p„ and pp are the number densities 
of the neutron and the proton and is the baryon density. In the present investigation, 
we study the effect of isospin asymmetry of the medium on the masses of the charmonium 
states J/ip,^{3686) and ^'(3770). 

The comparison of the trace of the energy momentum tensor arising from the trace 
anomaly of QCD with that of the present chiral model gives the relation of the dilaton field 
to the scalar gluon condensate. We have, in the limit of massless quarks 26| . 

= {^g;^^^^ ^ -(1 - d)x' (14) 

The parameter d originates from the second logarithmic term of equation (jTj). To get an 
insight into the value of the parameter d, we recall that the QCD /3 function at one loop 
level, for A^^^ colors and Nf fiavors is given by 

^.cM^~'-^(i-m^oi,^) (15, 



487r2 V lliVc 

In the above equation, the first term in the parentheses arises from the (antiscreening) self- 
interaction of the gluons and the second term, proportional to Nj, arises from the (screening) 
contribution of quark pairs. Equations f|T^ and f|T5l) suggest the value of to be 6/33 for 
three fiavors and three colors, and for the case of three colors and two fiavors, the value 
of d turns out to be 4/33, to be consistent with the one loop estimate of QCD /3 function. 
These values give the order of magnitude about which the parameter d can be taken 25| . 
since one cannot rely on the one-loop estimate for /3qcd(5')- In the present investigation of 
the in-medium properties of the charmonium states due to the medium modification of the 
dilaton field within chiral SU{3) model, we use the value of (i=0.064 [l^. This parameter, 
along with the other parameters corresponding to the scalar Lagrangian density, Co given 
by (!!]), are fitted so as to ensure extrema in the vacuum for the a, ( and x field equations, 
to reproduce the vacuum masses of the t] and rj' mesons, the mass of the a meson around 
500 MeV, and, pressure, p(po)=0, with po as the nuclear matter saturation density 

The trace of the energy-momentum tensor in QCD, using the one loop beta function 
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given by equation (|T5l) . for Nc=3 and Nf=3, is given as 

9a 



" 8 TT 

Using equations f|T4|) and f|T6|) . we can write 



= -o^^^M.G^^'^ (16) 



^G^G'^'^'^'^ = ^{1 - d)x' (17) 



We thus see from the equation (fTTI) that the scalar gluon condensate y^G^^^G^^'^'^j is propor- 
tional to the fourth power of the dilaton field, X) the chiral SU(3) model. As mentioned 
earlier, the in-medium masses of charmonium states are modified due to the gluon conden- 
sates. Therefore, we need to know the change in the gluon condensate with density and 
temperature of the asymmetric nuclear medium, which is calculated from the modification 
of the X field, by using equation (ITTl) . 

From the QCD sum rule calculations, the mass shift of the charmonium states in the 
medium is due to the gluon condensates [4, l8| . For heavy quark systems, there are two inde- 
pendent lowest dimension operators: the scalar gluon condensate ( {^G'^^^G^""''^ ) and the 
condensate of the twist 2 gluon operator ( {^^G^j^C'""^ ). These operators can be rewritten 
in terms of the color electric and color magnetic fields, {^E'^) and {^B"^). Additionally, 
since the Wilson coefficients for the operator {^B"^) vanishes in the non-relativistic limit, 
the only contribution from the gluon condensates is proportional to {^E"^), similar to the 
second order Stark effect. Hence, the mass shift of the charmonium states arises due to the 
change in the operator {^E"^) in the medium from its vacuum value jsj. In the leading 
order mass shift formula derived in the large charm mass limit the shift in the mass of 
the charmonium state is given as 8| 

^dipjk) ^2 k ( / (^s j^2\ /"^ 
k"^ /rric + e 

In the above, rrtc is the mass of the charm quark, taken as 1.95 GeV [8], is the vacuum 
mass of the charmonium state and e = 2mc — m^. ip{k) is the wave function of the charmo- 



nium state in the momentum space, normalized as / |^|V'(^)P = 1 |27|. At finite densities, 
in the linear density approximation, the change in the value of {^E'^), from its vacuum 
value, is given as 



1 0- 



and the mass shift in the charmonium states reduces to 

Am,(e) = -1 /^fc^l^r ,,/ (^E^) (20) 
' 9 J ^ dk 'P/mc + eXTT / n 2Mn 

In the above, (^^E'^'j^ is the expectation value of l^^E'^^ with respect to the nucleon. 

The expectation value of the scalar gluon condensate can be expressed in terms of the 
color electric field and the color magnetic field as 28 1 



'^Gl^G^''''^ = -2(^^{E^-B^)). (21) 



In the non-relativistic limit, as already mentioned, the contribution from the magnetic field 
vanishes and hence, we can write. 

Using equations ( 1T71) . ( fTSll and ( l22l) . we obtain the expression for the mass shift in the 
charmonium in the hot and dense nuclear medium, which arises from the change in the 
dilaton field in the present investigation, as 

In the above, x Xo are the values of the dilaton field in the nuclear medium and the 
vacuum respectively. 



In the present investigation, t 



le wave functions for the charmonium states are taken to 



be Gaussian and are given as [29 1 



= Normalization x Yi"\e, (j)){f3\^)^^exp-^^^'' L^^l^ (/^V^) (24) 

where = Mu/h characterizes the strength of the harmonic potential, M = mc/2 is the 
reduced mass of the charm quark and charm anti-quark system, and L^{z) is the associated 
Laguerre Polynomial. As in Ref. js], the oscillator constant /3 is determined from the mean 
squared radii (r^) as 0.46^ fm^, 0.96^ fm^ and 1 fm^ for the charmonium states J/'?/'(3097), 
'?/'(3686) and ■?/'(3770), respectively. This gives the value for the parameter /3 as 0.51 GeV, 
0.38 GeV and 0.37 GeV for J/^/'(3097), ?/'(3686 and ^(3770), assuming that these char- 
monium states are in the IS, 2S and ID states respectively. Knowing the wave functions 

9 



of the charmonium states and calculating the medium modification of the dilaton field in 
the hot nuclear matter, we obtain the mass shift of the charmonium states, J/V", '?/'(3686) 
and '?/'(3770) respectively. In the next section we shall present the results of the present 
investigation of these in-medium charmonium masses in hot asymmetric nuclear matter. 

III. RESULTS AND DISCUSSIONS 

In this section, we first investigate the effects of density, isospin-asymmetry and temper- 
ature of the nuclear medium on the dilaton field x using a chiral SU(3) model. From the 
medium modification of the x field, we shall then study the in-medium masses of charmo- 
nium states J/ip, ^/'(3686) and ^/'(3770) using equation fl23l) . 

The values of the parameters used in the present investigation, are : fco = 2.54, ki = 
1.35, k2 = —4.78, ^3 = —2.77, = —0.22 and d = 0.064, which are the parameters occurring 
in the scalar meson interactions defined in equation (jlj). The vacuum values of the scalar 
isoscalar fields, a and ( and the dilaton field x ^i-re —93.3 MeV, —106.6 MeV and 409.77 
MeV respectively. The values, ga-N = 10.6 and g(^iy = —0.47 are determined by fitting to 
the vacuum baryon masses. The other parameters fitted to the asymmetric nuclear matter 
saturation properties in the mean-field approximation are: g^jN = 13.3, gpp = 5.5, (74 = 79.7, 
gsp = 2.5, = 1024.5 MeV, m„ = 466.5 MeV and mg = 899.5 MeV. The nuclear matter 
saturation density used in the present investigation is 0.15 fm~'^. 

The variations of the dilaton field x with density, temperature and isospin asymmetry, 
within the chiral SU(3) model, are obtained by solving the coupled equations of motion of 
scalar fields given by equations ([9]), ( fTOl) . (fTTI) and ( fT2l) . In figure [H we show the variation of 
dilaton field x, with temperature, for both zero and finite baryon densities, and for selected 
values of the isospin asymmetry parameter, t] = 0, 0.1, 0.3 and 0.5. At zero baryon density, it 
is observed that the value of the dilaton field remains almost a constant upto a temperature 
of about 130 MeV above which it is seen to drop with increase in temperature. However, 
the drop in the dilaton field is seen to be very small. The value of the dilaton field is seen 
to change from 409.8 MeV at T=0 to about 409.7 MeV and 409.3 MeV at T=150 MeV and 
T=175 MeV respectively. The thermal distribution functions have an effect of increasing 
the scalar densities at zero baryon density, i.e., for /i*=0, as can be seen from the expression 
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FIG. 1: (Color online) The dilaton field x plotted as a function of the temperature, at given baryon 
densities, for different values of the isospin asymmetry parameter, rj. 

of the scalar densities, given by f|T3l) . This effect seems to be neghgible upto a temperature 
of about 130 MeV. This leads to a decrease in the magnitudes of scalar fields, a and (. This 
behaviour of the scalar fields is reflected in the value of which is solved from the coupled 
equations of motion of the scalar fields, given by equations ([9]), ( ITOl) . (|TT1) and (|T2l) . as a drop 
as we increase the temperature above a temperature of about 130 MeV. The scalar densities 
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attaining nonzero values at high temperatures, even at zero baryon density, indicates the 
presence of baryon-antibaryon pairs in the thermal bath and has already been observed in 
the literature 16|, l30[ • This leads to the baryon masses to be different from their vacuum 
masses above this temperature, arising from modifications of the scalar fields a and C,. 

For finite density situations, the behaviour of the x fi^ld with temperature is seen to be 
very different from the zero density case, as can be seen from the subplots (b),(c) and (d) 
of figure [H where the x field is plotted as a function of the temperature for densities po; 
2po and 4po respectively. At finite densities, one observes first a rise and then a decrease 
of the dilaton field with temperature. This is related to the fact that at finite densities, 
the magnitude of the a field (as well as of the C field) first show an increase and then 

n 

a drop with further increase of the temperature [1^ which is reflected in the behaviour 
of X field, since it is solved from the coupled equations of the scalar fields. The reason 
for the different behaviour of the scalar fields (a and Q at zero and finite densities can 
be understood in the following manner 16|. As has been already mentioned, the thermal 
distribution functions in ( |T3|) have an effect of increasing the scalar densities at zero baryon 
density, i.e., for /i*=0. However, at finite densities, i.e., for nonzero values of the effective 
chemical potential, pj*, for increasing temperature, there are contributions also from higher 
momenta, thereby, increasing the denominator of the integrand on the right hand side of the 
equation (fT3l) . This leads to a decrease in the scalar density. The competing effects of the 
thermal distribution functions and the contributions of the higher momenta states give rise 
to the observed effect of the scalar density and hence of the a and C, fields with temperature 
at finite baryon densities 16|]. This kind of behaviour of the scalar a field on temperature 
at finite densities has also been observed in the Walecka model by Li and Ko [31], which 
was reflected as an increase in the mass of the nucleon with temperature at finite densities 
in the mean field calculations. The effects of the behaviour of the scalar fields on the value 
of the X field, obtained from solving the coupled equations ([9]) to (fT2l) for the scalar fields, 
are shown in figure [H 

In figure [H it is observed that for a given value of isospin asymmetry parameter r], the 
dilaton field x decreases with increase in the density of the nuclear medium. The drop 
in the value of x with density is seen to be much larger as compared to its modification 
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with temperature at a given density. For isospin symmetric nuclear medium [rj = 0) at 
temperature T = 0, the reduction in the dilaton field x from its vacuum value (xo = 409.8 
MeV), is seen to be about 3 MeV at ps = Po and about 13 MeV, for pB = 4po- As we move 
from isospin symmetric medium, with r/ = 0, to isospin asymmetric medium, at temperature 
T = 0, and, for a given value of density, there is seen to be an increase in the value of the 
dilaton field x- However, the effect of isospin asymmetry of the medium on the value of the 
dilaton field is observed to be negligible upto about a density of nuclear matter saturation 
density, and is appreciable only at higher values of densities as can be seen in figure [H At 
nuclear matter saturation density, po, the value of dilaton field x changes from 406.4 MeV in 
symmetric nuclear medium (77 = 0) to 406.5 MeV in the isospin asymmetric nuclear medium 
{1] = 0.5). At a density of about 4po, the values of the dilaton field are modified to 396.7 
MeV and 398 MeV at = and 0.5, respectively. Thus the increase in the dilaton field x 
with isospin asymmetry of the medium is seen to be more at zero temperature as we move 
to higher densities. 

At a finite density, pb, and for given isospin asymmetry parameter rj, the dilaton field x 
is seen to first increase with temperature and above a particular value of the temperature, it 
is seen to decrease with further increase in temperature. At the nuclear saturation density 
Pb = Po and in isospin symmetric nuclear medium [r] = 0) the value of the dilaton field 
X increases upto a temperature of about T = 145 MeV, above with there is a drop in the 
dilaton field. For Pb=Po in the asymmetric nuclear matter with rj = 0.5, there is seen to 
be a rise in the value of x upto a temperature of about 120 MeV, above which it starts 
decreasing. As it has been already mentioned, at zero temperature and for a given value of 
density, the dilaton field x is found to increase with increase in the isospin asymmetry of the 
nuclear medium. But from figure [H it is observed that at high temperatures and for a given 
density, the value of the dilaton field x becomes higher in symmetric nuclear medium as 
compared to isospin asymmetric nuclear medium e.g. at nuclear saturation density Pb = Po 
and temperature T = 150 MeV the values of dilaton field x ^I's 407.3 MeV and 407 MeV at 
1] = and 0.5 respectively. At density pb = 4po, T = 150 MeV the values of dilaton field 
X are seen to be 399.1 MeV and 398.7 MeV for 77 = and 0.5 respectively. This observed 
behaviour of the x is related to the fact that at finite densities and for isospin asymmetric 
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FIG. 2: (Color online) The mass shift of i/ifj plotted as a function of the baryon density in 
units of nuclear matter saturation density at given temperatures, for different values of the isospin 
asymmetry parameter, rj. 

matter, there are contributions from the scalar isovector 5 field, whose magnitude is seen to 
decrease for higher temperatures for given densities, whereas 5 field has zero contribution 
for isospin symmetric matter. 

We shall now investigate how the behaviour of the dilaton field x the hot asymmetric 
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FIG. 3: (Color online) The mass shift of "0(3686) plotted as a function of the baryon density in 
units of nuclear matter saturation density at given temperatures, for different values of the isospin 
asymmetry parameter, r]. 

nuclear medium affects the in-medium masses of the charmonium states J/ip, ^'(3686) and 
?/'(3770). In figures [2], [3] and m we show the shifts of the masses of charmonium states J/ip, 
'0(3686) and '0(3770) from their vacuum values, as functions of the baryon density for given 
values of temperature T, and for different values of the isospin asymmetry parameter, t]. 
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FIG. 4: (Color online) The mass shift of "0(3770) plotted as a function of the baryon density in 
units of nuclear matter saturation density at given temperatures, for different values of the isospin 
asymmetry parameter, r]. 

We have shown the results for the values of the temperature, T = 0, 50, 100 and 150 MeV. 
At the nuclear matter saturation density, = Po at temperature T = 0, the mass-shift for 
J/ip meson is observed to be about —8.6 MeV in the isospin symmetric nuclear medium 
[r] = 0) and in the asymmetric nuclear medium, with isospin asymmetry parameter rj = 0.5, 
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it is seen to be —8.4 MeV. At pB = 4po, temperature T = 0, the mass-shift for J/V' meson 
is —32.2 MeV in the isospin symmetric nuclear medium (r] = 0) and in isospin asymmetric 
nuclear medium [t] = 0.5), it changes to —29.2 MeV. The increase in the magnitude of the 
mass-shift, with density pb, is because of the larger drop in the dilaton field x higher 
densities. However, with increase in the isospin asymmetry of the medium the magnitude 
of the mass-shift decreases because the drop in the dilaton field x is less at a higher value 
of the isospin asymmetry parameter rj. At the nuclear matter saturation density pb = Po, 
and for temperature T = 0, the mass-shift for ■?/'(3686) is observed to be about —117 and 
— 114 MeV for values of the 77 = and 0.5 respectively, and for ■?/'(3770), the values of 
the mass-shift are seen to be about —155 MeV and —150 MeV respectively. At pb — 4po 
and zero temperature, the values of the mass-shift for ■0(3686) are modified to —436 MeV 
and —396 MeV for 77 = and 0.5 respectively, and, for ■0(3770), the drop in the masses 
are about —577 MeV and —523 MeV respectively. As mentioned earher, the drop in the 
dilaton field, x, at finite temperature is less than at zero temperature and this behaviour 
is reflected in the smaller mass-shifts of the charmonium states at finite temperatures as 
compared to the zero temperature case. At nuclear matter saturation density pb — Po, 
and for temperature T — 100 MeV, the values of the the mass-shift for the J/t/j meson 
are observed to be about —6.77 MeV and —6.81 MeV for isospin symmetric (77 = 0) and 
isospin asymmetric (77 — 0.5) nuclear medium respectively. At baryon density pb — 4po, 
temperature T = 100 MeV, the mass-shift for J/ip is observed to be —28.4 MeV and —27.2 
MeV for isospin symmetric {rj = 0) and isospin asymmetric [r] = 0.5) nuclear medium 
respectively. For the excited charmonium states ^0(3686) and ■0(3770), the mass-shifts at 
nuclear matter saturation density Pb — Po and temperature T = 100 MeV, are observed to 
be —91.8 MeV and —121.4 MeV respectively, for isospin symmetric nuclear matter (77 = 0) 
and —92.4 MeV and —122 MeV for the isospin asymmetric nuclear medium with 77 = 0.5. 
For baryon density ps = 4po, and temperature T = 100 MeV, mass-shift for the charmonium 
states -0(3686) and 0(3770) are modified to about —386 MeV and —510 MeV respectively, 
for isospin symmetric (77 = 0) and —369 MeV and —488 MeV for isospin asymmetric nuclear 
medium with 77 = 0.5. 

For temperature T — 150 MeV and at the nuclear matter saturation density pb — Po, 
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the mass-shifts for the charmonium states J/ilj,tp{3686) and iIj{3770) are seen to be —6.25, 
—85 and —112 MeV respectively in the isospin symmetric nuclear medium [r] = 0). These 
values are modified to —7.2, —98 and —129 MeV respectively, in the isospin asymmetric 
nuclear medium with rj = 0.5. At a baryon density pB = 4po, the values of the mass-shift 
for J/ip,ilj{3686) and ^/^(STTO) are observed to be —26.4, —358 and —473 MeV in isospin 
symmetric nuclear medium [rj = 0) and in isospin asymmetric nuclear medium with rj = 0.5, 
these values are modified to —27.6, —375 and —494 MeV respectively. Note that at high 
temperatures e.g at T = 150 MeV the mass-shift in isospin asymmetric nuclear medium 
[rj = 0.5) is more as compared to isospin symmetric nuclear medium [rj = 0). This is 
opposite to what is observed for the zero temperature case. The reason is that at high 
temperatures the dilaton field x has more drop in the isospin asymmetric nuclear medium 
[r] = 0.5) as compared to the isospin symmetric nuclear medium (77 = 0), due to the 
contributions from the 6 field for nonzero rj, which is observed to decrease in its magnitude 
at high temperatures. 

The values of the mass-shift for the charmonium states obtained within the present in- 
vestigation, at nuclear saturation density po and temperature T = 0, are in good agreement 
with the the mass shifts of J/^, ^(3686) and ^(3770) as -8, -100 and -140 MeV re- 
spectively, at the nuclear matter saturation density, computed in Ref. js] from the second 
order QCD Stark effect, with the gluon condensate in the nuclear medium computed in the 
linear density approximation. The mass-shift for J/ ip has also been studied with the QCD 



sum rules in jsl and the value at nuclear saturation density was observed to be about —7 
MeV. In Ref. 3| the operator product expansion was carried out upto dimension six and 
the mass shift for J / ip was calculated to be —4 MeV at nuclear matter saturation density po 
and at zero teraperature. The effect of temperature on the J/ip in the deconfinement phase 



was studied in 
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32[. In these investigations, it was reported that J/ip mass is essentially 



constant in a wide range of temperatures and above a particular value of the temperature, 
T, there was observed to be a sharp change in the mass of J/ip in the deconfined phase 



33| . In the present work, we have studied the effects of temperature, density and isospin 



asymmetry, on the mass modifications of the charmonium states {J/ip, '?/'(3686) and '?/'(3770)) 
in the confined hadronic phase, arising due to modifications of a scalar dilaton field which 



18 



simulates the gluon condensates of QCD, within a chiral SU(3) model. The effect of tem- 
perature was found to be small for the charmonium states J/'?/' (3097), ?/'(3686) and iIj{3770), 
whereas, the masses of charmonium states observed to vary considerably with density, in 
the present investigation. 

In summary, in the present work, we have investigated the effects of density, temperature 
and isospin asymmetry of the nuclear medium on the masses of the charmonium states 
J/V', ^'(3686) and '?/^(3770), arising due to modification of the scalar dilaton field, X; which 
simulates the gluon condensates of QCD, within the chiral SU(3) model and second order 
QCD Stark effect. The change in the mass of J/tp with density is observed to be small at 
nuclear matter saturation density and is in agreement with the QCD sum rule calculations. 
There is seen to be an appreciable drop in the in-medium masses of excited charmonium 
states '?/'(3686) and '?/'(3770) with the density. At the nuclear matter saturation density, the 
mass shifts of these states are similar to the values obtained using the QCD second order 
Stark effect with the modifications of the gluon condensates computed in the linear density 
approximation For a given value of density and temperature, the effect of the isospin 
asymmetry of the medium on the in-medium masses of the charmonium states is found to be 
small. This is due to the fact that the magnitude of the S field remains small as compared to 
the a and ( fields. At finite densities, the effect of the temperature on the charmonium states 
is found to decrease the values of mass-shift upto a particular temperature, above which the 
mass shift is seen to rise. This is because of an initial increase in the dilaton field x and then a 
drop with further increase in the temperature, at a given baryon density, arising from solving 
the coupled equations for the scalar fields. This is related to the fact that the scalar densities 
of the nucleons initially drop and then rise with the temperature at finite values of the baryon 
densities. The mass drop of the excited charmonium states (\l/(3686) and ^(3770)) are large 
enough to be seen in the dilepton spectra emitted from their decays in experiments involving 
p-A annihilation in the future facility at GSI, provided these states decay inside the nucleus. 
The life time of the J/\l/ has been shown to be almost constant in the nuclear medium, 
whereas for these excited charmonium states, the lifetimes are shown to reduce to less than 
5 fm/c, due to appreciable increase in their decay widths 29(|. Hence a significant fraction of 
the produced excited charmonium states in these experiments are expected to decay inside 
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the nucleus |34|- The in-medium properties of the excited charmonium states ?/'(3686) and 
iIj{3770) can be studied in the dilepton spectra in p-A experiments in the future facihty of 
the FAIR, GSI. The mass shifts of the charmonium states in the hot nuclear medium seem 
to be appreciable at high densities as compared to the temperature effects on these masses, 
and these should show in observables like the production of these charmonium states, as well 
as of the open charmed mesons in the compressed baryonic matter (CBM) experiment at the 
future facility at GSI, where baryonic matter at high densities and moderate temperatures 
will be produced. 
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